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SUMMARY 

In this report, we address tfie problem of pilot signal parameter 
optimization and the related problem of power transponder performance 
analysis for the SPS reference phase control system. Signal and 
interference models are established to enable specifications of the 
RF front end filters including both the notch filter and the antenna 
frequency response. A simulation program package is developed to be 
included in SOLARSIM to perform tradeoffs of system parameters based 
on minimizing the phase error for the pilot phase extraction. 

An analytical model that characterizes the overall power transponder 
operation is developed. From this model, the effects of different phase 
noise disturbance sources that contribute to phase variations at the 
output of the power transponders can be studied and quantified’. 

The important findings on the transponder design parameters 
and results can be summarized as follows: 

©EIRP = 93.3 dBW 

-oPN Chip Rate ^ 10 Mcps 

©RF filter 3 dB, cutoff frequency 20 MHz 

•Notch filter 3 dB cutoff frequency % 1 MHz 

©Notch filter dc attenuation ^ 60 dB 

•©PN Code period > 1 msec 

©Costas loop phase jitter £0.1 deg for 10 Hz loop bandwidth 
©PN Code loop jitter is negligible 

/ 

©Channel Doppler is negligible 

©Klystron phase control loop bandwidth > 10 KHz 
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In conclusion our results indicate that it is feasible to hold the 
antenna array phase error to less than one degree per power module 
for the type of disturbances modeled in this report. However, there 
are irreducible error sources that are not considered, herein and their 
effects remain* to be seen. They include: 

©Reference phase distribution errors 
©Differential delays in the RF path 
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1.0 INTRODUCTION 

This report serves to document results from the Pilot Signal 
Parameter Optimization/Analysis (Task I) and the Power Transponder 
Analysis/Modeling (Task II) of the SPS Antenna Phase Control System 
Hardware Simulation Study (Phase III). It can be divided into two 
parts: Sections 2.0-11.0 deal with Task I and Sections 12.0 - 15.0 

are devoted to Task II. 

The key technical issues to be addressed under Task I can be 
summarized by Fig. 1.1. The interferences are generated by the following 
mechanisms: (1) self jamming due to the power beam leakage from the 

diplexer/circulator; (2) mutual coupling from adjacent transponders 
and (3) thermal noise. The signal and interference spectrum at the 
input to the SPS* transponder is depicted in Fig. 1.1(b). in general, 
the combined phase noise inference from the power beams consists of a 
coherent and a noncoherent term. Depending on the mechanization of 
the antenna structure and diplexer/circulator characteristics, these 
terms are associated with gains K-j and K^. Note that the phase noise 
interferences are concentrated around the carrier frequency (2450 MHz). 
The uplink pilot signal on the other hand has no power around this 
frequency. Its power spectrum peaks at f % 0.75 R c , with a value 
proportional to the product of the received power (P^) and the PN 
chip rate (R c ), and inversely proportional to the PN code length (M). 

The parameters R c and M are also related to the processing gain of the 
PN spread signal and determine its interference suppression 
capability. Notice also that about 14% of the pilot signal power 
lies between [f| < 2R c . The RF filter characteristic is mainly 
determined by the waveguide antennas, which have bandwidths ranging 
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Figure 1.1. Signal and Noise Spectrum into SPS Transponder. 
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from approximately 15 to 45 MHz depending on the array area. Our 
goal is to optimally select (1) the pilot signal so that it passes through 
the RF filter with negligible distortions, and (2) a practical notch 
filter that rejects most of the phase noise interferences. When this 
is done, one can De assured that the reconstructed pilot signal phase 
after the sync loops bears a tolerable error for the retrodi recti ve 
scheme. 

In this report, we have characterized the interference model and 
requirements for the optimization of the first IF filter for the pilot 
signal communication system. We have characterized analytically the 
power spectral density of the pilot signal as well as various sources 
of interference. From this information, we are able to optimize 
the first IF filter in the reference SPS transponder for interference 
rejection. We have also formulated a mathematical framework which 
serves as a basis where different tradeoffs can be made 'in terms of 
system parameters such as pilot signal transmitter EIRP, PN code 
requirements and chip rates. Based upon this mathematical model, a 
computer program is developed to be included in the SOLARSIM package to 
perform tradeoffs of pertinent design parameters of the receiver portion 
of the SPS transponder; the phase error of the pilot phase tracking 
(Costas) loop is used as the performance measure. As a result of the 
SOLARSIM data, the following can be specified on a preliminary basis: 

(1) chip rate, (2) notch filter response, (3) code period, (4) uplink 
EIRP and (5) Costas loop bandwidth. 

Under Task II, we have developed analytical models for the SPS 
transponder tracking loop system that includes: (1) the PN despreader 
loop, (2) the pilot phase tracking (Costas) loop and (3) the PA phase 
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control loop. The phase reference receiver that feeds the reference phase 
distribution system is also modeled. Various sources of potential 
phase noise interferences are identified and their effects on the 
performance of the individual loops are modeled. In particular, a 
model of the phase noise profile of the klystron amplifier based 
on a specific tube measurement is introduced. Important implications 
on the PA control loop design are also addressed. 

An analytical model for evaluating the overall performance of the 
SPS transponder is given. The phase fluctuation at the output of the 
transponder is shown to be directly related to the various noise 
processes through the closed-loop transfer functions of the tracking 
loops. THese noise processes are either generated externally to the 
transponder circuitry such as ionospheric distrubances, transmit 
frequency instability, or externally such as receiver thermal noise, 
power beam interferences, data distortions, VCO/mixer phase noise 
and the phase variations introduced by the reference distribution tree. 

A detailed computer simulation is deemed necessary to quantitatively 
investigate the interplays between the .elements of the transponder. 
Based on our preliminary investigations, some recommendations on the 
transponder design are also made. 
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2.0 THE PROBLEMS AND ASSUMPTIONS ASSOCIATED WITH SPS PILOT SIGNAL 
DESIGN OPTIMIZATION 

We wish to design the pilot signal communication system to operate 
reliably in the face of several types of interference: 

(1) The downlink power beam signal. 

(2) Noise in the spacetenna receivers. 

(3) Unintentional or intelligent RFI. 

(4) Intelligent beam stealing signals. 

Figure 2.1 illustrates the R.F. signal scenario under consideration. The 
spacetenna is -composed of TOT, 552 variable size rectangular arrays ’which range 
from 1.73mxl.73m to 5.2mx5.2m, operating at a frequency of 2450_MHz. 

The power beam signal contributed by the i— subarray is denoted 
by s^t) for i=l,2,...,N, measured at the output port of the i— - 
transmitter. The signal at the output of the j — - subarray's RF 
receiver is given by r. (t) and possesses components due to the 

J 

power beam signals s'.(t), i=l , ,N, the receiver noise n.(t), the 

1 J 

RFI s R pj(t), and a possible beam-stealing signal s BS (t), in addition 

to the desired pilot signal s p (t). If we knew all of these input 

signals to the scenario, we could develop a representation for r.(t) 

3 

as indicated in Figure2.2 ,once we have the following additional 
information : 

(a) The system function H. .(f) describing the coupling of the 

* 

XL XL 

i— power beam signal to the j— - receiver's RF output signal 
r. (t). This includes the effects of all waveguides, 

v 

circulators, RF receiver filters, antennas, etc. This 
must be known for all values of i. 

(b) The system function H RF (f ) from the j— -subarray 's antenna 

X L 

terminal to the j— receiver's RF output signal r-(t). This 

\J 
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Figure 2.1. RF Signal Scenario.. 
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also includes the effects of waveguides, circulator, RF 
filters, etc. 

(c) Atmospheric and ionospheric channel models for the paths 

from the pilot signal transmitter, the RFI source, and 

the beam stealer to the antenna terminal of the j— subarray. 

< 

(d) A “threat model" indicating the capabilities and level of 
sophistication of the beam stealing processor. 

Knowledge of the above quantities would make the analyses to 
follow more precise. 

The pilot signal baseline design [1] is given as NRZ/BPSK/ 
Bi- 4 -DS/CDM, i.e., 

Sp(t) = Re{d(u,t)c(u,t)a(u,t)e Ja)ct } 


where u is a sample space parameter indicating randomness, d(u,t) ts 
an NRZ binary data signal BPSK modulated on a Manchester-coded 
direct-sequence spread spectrum signal, and a(u,t) represents 
(unwanted) amplitude and phase modulation on the pilot transmitter. 
We assume that the following approximations are valid: 

(a) The atmospheric-ionospheric channel can be modelled as: 


s D (t) 


MOD. bH 


H CH ( f ) 


(SIGNAL AT PILOT f 
TRANS. OUTPUT) b(u,t) 


Sp(t) 


(SIGNAL AT SUBARRAY 
APERTURE) 


3« r t‘ 

Sp(t) = j h CH (t-t')Re{d(u,t , )c(.u,t , )a(u,t')b(u,t , )e }dt' 

—Co 

(b) The filtering H CH (f) imposed by the channel is 1 wide 

compared to- the receiver RF filtering H^p(f). Assuming 
that attenuation functions which are uniform across the 
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bandwidth of H R p(f) are included in the channel fading 
function b(u,t), we approximate h CH (t) 


h CH (t) = $ (t-r(u)} 

by a pure random delay. Therefore 

j(o_t 

Sp(t) = Re{m * p Ct-x (u) )e c } 
m p(t) = a(u,t)b(u ,t)c(u,t)d(u,t) 


(c) The RF and IF sections of the subarray receivers are 

modelled in Figure 2.3. Under weak narrowband assumptions 
we can reduce the RF and IF subsystem signal processing 
model to an equivalent complex baseband signal processing 
model. These techniques are discussed in Appendices 1 
and 2. The resultant model is shown in Figure 2.4. If s(t) 
is the input signal in Figure 2.3, then m (t) is the 

A 

equivalent input signal in Figure 2.4 where 

3a t 

'->! 'x(t) = Re{m x (t)e c } 

The output signal y(t) and its baseband equivalent m (t) 

* */ 

ace related by 

jw-rpt 

y(t) = Re{m (t)e } 


where the IF frequency is 490 MHz which is one-fifth 
of the carrier frequency. 

(d) We further assume that the factor a(u ,t)b(u ,t}d(u,t) 

of the received pilot signal's modulation is quite narrow 
band relative to the passband of the equivalent RF system 
function G RF (f). The product filtering approximation of 
Appendix 3 can then' be applied to further simplify the effect 
of the RF filter on the received pilot signal. 




C inborn 



64 » 



am 


OUTPUT 
TO PN 

SYNC SYSTEM 


INPUT 

INCIDENT 
ELECTRIC 
FIELD IN 
SUBARRAY 
APERTURE 


RF SYSTEM 
FUNCTION 


H RF (f) 


IF 

MIXER 


FIRST IF 
FILTER 


H-j (f ) 


cos(.8co c t+0(u,t)) 

FROM CARRIER SYNC 
SYSTEM 


SECOND IF 
FILTER 


H P (f) 


OUTPUT 

TO CARRIER~SYNC S 
AND DATA DEMOD 
SYSTEMS 


C R (u,t-x(u)) 

FROM PN SYNC 
SYSTEM 


Figure 2.3. RF and IF Subsystem Signal Processing. 







8S20 S& 



TO SYNC 



G RF (f) « U(f + f c )H RF (f + f c ) 

6-,(f) = U ( f +f j F ) H -j (f+f IJr ) 
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Figure 2. 4. Complex Baseband Equivalent Circuit for the RF and IF Processor. 
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3.0 POWER SPECTRAL DENSITY COMPUTATIONS 


We will now sketch the computation of the power spectral density 
at the input to the second IF filter. This is useful information 

for determining carrier tracking loop and 'data detector performance. 
The PN tracking loop contains a similar signal and these computa- 
tions can be adapted to those design needs as well. We assume 
that superposition is valid and analyze each signal source's effect 
separately. We feel that this is a very good assumption despite the 
number of tracking loops which feed very slowly varying functions 
of the IF signal back into the system. These loops are "broken" 
in our computation. 

3.1 The Pilot Signal 

The Bi~ 4 > DS spread spectrum (SS) modulation c(u,t) on the pilot 
signal is of the form 

c(u,t) = a i P(t-iT c -6 (u)) 

i 

where a., is a known +1 binary sequence, S (u) is a random delay, T c is 
the chip time, and P(t) is the basic Bi-<}> waveform: 



The chip time T c and the sequence La..} are design variables. 

Pure power spectral analysis cannot be used to analyze the 
effect of the RF filter on the pilot signal. This is based on the 
fact that the phase of the SS code on the pilot signal relative to 
the phase of the code reference signal from the PN tracking loop 
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determines the bandwidth occupancy of this pilot signal component of 

the input to the second IF filter. 

When Sp ( t ) passes through the RF and first IF filters, we 

assume that the main effect of these operations is the reshaping 

, of P(t). Hence the output of the first IF filter is 

where 

= atu.tjbtu*): 0 ) (u,t)d(u,t)e Je(u ’ t+r(u)) 

c (1) (u,t) = 2 a.Qtt-IT -s(u)) . 
i ^ 

and the new pulse shape Q(t) is given by 

Q(t) = P(t)*g RF (t)*g 1 (t) . 

Here g R p(t) and g-j(t) are the impulse responses of the first two 
filters in the complex baseband equivalent circuit (Figure 4), and 
* denotes convolution. 

After PN code reference multiplication, the baseband equivalent 
input to the second IF filter is: 

m^(u,t) = m^(t-T(u))c R (u,t-x(u)) 

= A(u,t)c^^(u,t-T(u))c R (u,t-T(u)) 

where A(u,t) now represents the low bandwidth signal: 

A(u,t) = a(u,t-T(u)}b(.u,t- T (u))d(u,t-T(u))e :i0(lJ ’ t) 

(We in fact assumed earlier that A(u,t) has .bandwidth « 1/T .) 

c 

Now c R (u,t) is a replica of c(u,t) except that the transmitter 
clock error s (u) is not known. 
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. 1 = 53 -a.Pft-i^) 

1 

The product of SS signals appearing in m^(u,t) is* then 
C^(u,t) = C^(u ,t-r(u))c R (U s t-r(u)) 

i 

x 53 ajPft-jTc-rfu)) 

3 

- 53 f)(u,t) 
k 

« 

where 

c[ 2 ){u,t) = 23 a j a j+k^( t- ( J+k ) T c~ 6 ( u )" T ( u ))P( t “^ T c"^( u )) 

3 

In the ideal processor from the pilot signals point of view, the 
RF filter does not distort the pulse, i.e. Q(t) = P(t), the PN sync 
subsystem adjusts t(u) = t(u)+6(u), and hence we would have 

c^(u,t) = Cq 2 ^ (u , t) = 1 (Ideal Pilot Signal Reception) 

Due to RF filter distortion and PM tracking error, Cg 2 ^(u,t)‘ will 
be a periodic waveform with period T . One period of the waveform 
{between -T 72 + t(u) and T 72 + T (u))is given by P(t- x (u))Q(t-s (u) 
-t(u)). Lets expand this waveform in a Fourier series to determine 
its harmonic content: 


Cq 2 ^ (u jt ) = ^ P n (u)exp(j2irnt/T c ) 


tnx^om 




T c /2+;(u) 


P(t- T (u))QCt-6(u)- T (u))expC-j27rnt/T c )dt 


-T c /2+ x (u) 


exp 


( j2irnx(u)^ 

V. T c J 


[ T c /Z 

i-T c /2 


P(t)Q(t-e(u))exp(-j2TTnt/T c )dt 


and the PN apparent tracking offset is given by 


e(u)' = x(u) + 5 (u) - xCu) 


We call this the apparent tracking error for two reasons: 

(1) Ideally we would like the PN tracking loop to drive 
e(u) to the value which maximizes |pg(u)j to get the 
maximum energy into the narrow 2 nu IF filter and later 
detection filters. This maximum may not occur at e(u)=0. 

(2) The fact that P(t)^Q(t) and Q(t) is not generally symmetric 
in any way implies that the "S-curve 11 of the PN tracking 
loop will not generally be symmetric or have a zero at 
e=0. 

As indicated in (1) above, we have identified a portion of the pilot 

n j 

signal which corresponds to the desired 2— IF filter input, namely 
the signal whose complex baseband representation is given by 

mp 3 ^(u,t) = c^(u,t)A(u,t) 

* P 0 (u)A(u,t) 

The last approximation is quite accurate as far as 2— IF filter 

(21 

computations go, since the other spectral components of Cq ; (u,t) 
shift the spectrum of A(u,t) away from the filter passband by integer 
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multiples of 1 /T c - 

(21 (21 

The remaining components cjj. J [u,t), k^O, of '(u,t) 
correspond to signals modulated at the chip rate by the sequences 
^ a j a j+k} ’ anc * * lence are i 1 1 spread spectrum signals. We view 
these as distortion noise and define their baseband equivalent 
signal as 

n‘ 3 W) & A(u.t) E c£ 2 >(u,t) 

p k^O K 

The sum signal has period MT c where M is the period of the sequence 
{a^} and can be analyzed spectrally once the design parameters are 
selected. For now we assume that the power in n^(u,t) which is 

r 

in the passband of the IF filter is lower than that of m^(u,t) by 

at least the processing gain of the SS system. Hence is is not a 

major factor in the preliminary design and will be ignored. 

At this point power spectral density calculations are quite 

easily done for the equivalent baseband pilot signal m^(t). Let 

P 

us generally define 


R x (t) = <E{x(u ,t+T )x*(u ,t) }> 


and 


S v (f) - 


■x- ■ / X 


R (t ) exp ( - j 2nf t ) dr 


We assume that the data modulation d(u,t) is independent of B(u,t) 
where 

B(u,t) = a{u,t)b(u,t)e j0 ( u,t+T ( u ^ 


But e(u,t) cannot be considered independent of a(u,t) and b{u,t) 
whose phases it is designed to track. Continuing then, 
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nip 3 ^(u,t) = P 0 (u)B(u,t)d(u,t) 
and 

S ( 3 )( f ) = |P 0 ( U ) | 2 >S B ( f )* s d ( f ) 

m p 

3.2 Independent Interference 

One of the major advantages of SS systems is their relative 
immunity to wide-sense stationary interference. There are several 
sources of interference which are independent of the SS code 
reference signal c^(u,t): 

1. Receiver noise nj(u,t). 

2. Unintentional RFI (we exclude multipath and intelligent 
jamming from this category). 

3. Downlink signals s-(t), j=l,2,...,N. 

J 

The theory of wide-sense stationary independent interference rejection 
is covered in [2]. 

In general if a narrowband signal s(u,t) is represented by a 

complex baseband equivalent signal m s (u,t), i.e., 

Jut 

s(u,t) = Re{m $ (u,t)e c } 

Then the power spectral density of m s (u,t) can be determined from 
the corresponding density of s{u,t) by 

* 4U(f + f 0 )S s (f + f 0 ) 

or, reversing the relation, 

S s (f > “ k (f -' f o> + V ( - f - f o ) 3 

s s 

We assume that we can determine the following power spectral 
densities : 
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1. S (f) = N n /2 = equivalent noise power spectral density* 

n j U 

of the receiver reference to the antenna terminals. 

2. S RF i(f) = the RFI power spectral density measured at the 
subarray antenna terminals. We assume that the effects 
of the RFI channel have been taken into account already. 

3. S (f) - the power spectral density of the output of the 

S 3 

j-th subarray transmitter. 

4. c ( f ) = the power cross-spectral density of the outputs 

s i i> k 

J K th th 

of the j — and k — subarray transmitters. This is the 
Fourier transform of 

R s.s = < ^ s jCu;t+T)s^Cu,t)}> 

j k 

We expect these to definitely be nonzero since all subarray 
transmitters are transponding the same pilot tone with, 
basically the same modulation on it. 

We further assume that for purposes of this spectral computation 
we can set the carrier tracking phase e(u,t) to a constant e(u). 

Now it is a simple matter to compute the power spectral density 
of the interference at the input to the second IF filter, due to these 
signals. The resultant equivalent baseband power spectral density is: 


S m (f) ' s c (f >* 

I C R 


S,(f) 


N N 

+ Z S s 

i=l k»l 


(2N 0 +S V i (f))|G RF (f) 


Vk (f >V f)G V f > 


where 


\FI <f) = 4U(f+ VW f+ V 

* W +f 0 > S s. S ^ 

i K ^ IK 
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and 

V, (f) = ' 4u l f4 V s s 1 < f4 V> 

The function S (f) represents the power spectral density of the 
C R 

PN code reference signal generated in the receiver. 

3.3 Intelligent Jamming 

In this situation the jammer attempts to emulate the pilot tone, 
either to add narrowband- noise modulation and jam the data signal, or 
to heam steal. To be effective in either case, the jammer must 
initially synchronize his SS code modulation with the pilot tones 
SS code modulation so that they arrive nearly in phase at the 
receiver. He must also deliver a signal power level at least 
comparable to the pilot signal’s level at the spacetennas subarray 
terminals. The basic power density analysis is similar to that used 
for the pilot signal earlier. 


4. DESIGNING THE FIRST IF FILTER 


In this portion of the design it seems reasonable to assume that 
the unintentional RFI and the intelligent jamming signals are not 
present. The SS processing (code multiplier and second IF filter) 
is the key to rejecting these interference signals. The major 
purpose of the first IF filter is to reject as much of the extremely 
strong downlink signal -as is possible. 

One plausible criterion for the design of H^(f) or its baseband 
equivalent G-j(f) is to maximize the signal -to -interference density 
ratio SIDR in the center of the second IF filter passband. Thus we 
consider 


' cyCinC^om. 
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f{jP 0 ^! > S B (f) * S d (f) 

SI DR - *— ~n " ~ 

S c Cf)*[|6,(f)| SjCf)] f=0 
R , 


where 


T c /2 

P 0 (u) = .p f P(t)Q(t-e(u))dt 

C J - T c /2 

which by Parsevals theorem is 

P°C) ' ^ f |p(f)| 2 S RF (f) e ' 3 “ e(u) G, <f ) df 

— co 

with p(f) being the Fourier transform of P(t). Also the quantity 
S r (f)* in the denominator is given by 

V f ) * 2N ol fi RF' f )| 2 + i: E W f)G 1i (f)e & (f) 

1=1 k=l 1 K 

Another plausible criterion for the design of the first IF filter 
is the maximization of the signal-to-interference power ratio in the 
passband of the second IF filter. 


SIPR = 


CO- 

t{|P 0 (u)| 2 }/ |G 2 (f)| 2 [S B (f)*S d (f)]df 

-00. 

r"|G 2 (f)| 2 (S CR [f)*[|G 1 (f)| 2 S I (f)])df 


The problem with both of these performance measures is their dependence 

on parameters which are not decided at this point. 

Before continuing we must get an idea of the shape of S (f). Its 

C R 

spectral density is computed in Appendix 4. Below is a sketch of the 

envelope of the line spectrum S (f) assuming that the spectrum S (n) 

c R - a 

of the sequence {a^} is flat (see Appendix 4 for definitions.) 


cjClnC^oi 


inx^om 


s r (f) 


Envelope of Line Powers 
■in S c (f). Line Spacing=l/MT C 

M = P§riod of la ). 

n 


When the line spacing 1/MT in S (f) is much less than the 

c c R 

width of Sjff), and assume S & (f) is flat, we can approximate 


S* (f)*[jG (f)|^S T (f)] 

C R 1 1 f=0 


constant xj |p(f)G-j (f)‘pSj(f)df 


and hence 


SIDR % constant x 


;|J" |p(f)| 2 6 RF (f)e' 3 “ e(u) e 1 tf)df 

f |p(f)G n (-F) l 2 Sj(f)df 

•f —co 


In many SS systems, the line spacing 1/MT is also the data 

V Si 

nd 

bandwidth and the approximate bandwidth of the 2— IF filter. 
In this case we make the approximation 


f” |G 2 (0| 2 S c (f-a)df 

m 


constant x jp(a) |‘ 


cJLinC^oi 
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where |p(f)r is the envelope of the line spectrum S (f), and 

R 

again we assumed S fl (n) is flat. This can be used to significantly 
simplify the denominator of SIPR, and 


SIPR ' u constant x 


{iT |P(f)| 2 G RF (f > e " J “ e(U)G 1 <f)df 2 | 

f (-F) | z Sj(f)df 

♦'-CO 


So under certain circumstances both SIPR and SIDR have approximately 
the same dependence on G-j(f) shown explicitly above. 

Lets work on the numerator expression to evaluate the expected 
value and determine its relation to the denominator. 


Numerator = 


jj |P(f-j)| ®Rp(f 1 )^*(f i ) J P (f 2 ^ I G RF^2^ G 1^2^ 


j(a)-,-0) ? )elu) 

E{e 1 >df-jdf 2 


To continue precisely much depends on the characteristic function 
of the tracking error. The integrals factor only if we assume 
e(u) % £{e(u)} = m e , i.e. e(u) is very close to its mean value 
with probability approaching one. Lets assume this is true and 


see what happens: 


Numerator = 


-julll 1 2 


lP(f)re RF (f)G 1 {f)e e df 


Applying Schwartz's inequality then gives 


Numerator 


< j |g 1 (f ) P (f )v^rm 2 df 

— CO 

'• » i — joim } 2 

‘ X j Aj(f)' G RF< f >P*^ e 


P- f 7 
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or 


Numerator f “ 
Denominator — j 

co 


lG RF (f)p(f)| 2 

K - y - m — df * SIDR or SIPR 

max max 


The maximum on the right is independent of G^(f) and is achieved 
when 


Jwni 


We 


%(f)e 
S,( f ) as jprj 

must now try to design a decent approximation to G-j(f), even 


though it may be difficult to approximate m and determine S T (f) 

6 1 

precisely. 

5.0 APPROXIMATIONS TO THE INTERFERENCE SPECTRUH S ^f) 

The major problem encountered in attempting the design of H-j(f) 
(or the equivalent G^(f}) is the evaluation of Sj(f). Recall that 


s j( f ) - 2N 0 |6 RF (f)|^ + 


S m (f)G.,.(f)G*.(f) 
ITI i m k lj 


i k 

represents the equivalent baseband power spectral density of the 
interference at the RF output of spacetenna receiver j. Hopefully 
hardware experience will provide us with estimates of N^/ 2 and G R p(f}. 
The complicated double sum however requires some simplification. 

We shall write 1 the following expression for the signal transmitted 
by the k— antenna subarray. 

ju t 

s k (t) = Re{m k (u,t)e } 


■where 


j[<j>{u,t-d k )+o) c d k ] 


m k (u,t) = A k (u,t)e 


d k = The differential delay between the pilot tone arrival 
at subarray k and the pilot tone arrival at subarray j. 


'cjCinC^orri 
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th 

<!>(u,t) = The phase function appearing on the output of the j— 
subarray transmitter in the absence of transmitter 
noise and receiver noise. The time variations in <f>(u,t) 
are due to the modulation factor a(u,t)b(u,t) on the 
pilot tone. 

“ Differential carrier phase shift between the pilot tone 
at subarray k and the pilot tone at subarray j. 

A^Cu , t ) = Effects of transmitter and receiver noise generated in 
subarray k. 

Ideally A^(u,t) = a^ for all k, and m k (u,t) would be the same as 
m- (u,t) except for a time delay and a constant complex phase factor. 

J * 

We now make the following assumptions: 

(a) {Aj,(u,t)> is a set of independent, identically-distributed 
wide-sense stationary random processes each with mean 

and covariance K^(t). 

E{A k (u,t)> = .A, Vk,t . 
CovWjJU'st-jM^u,^)} = K A (t r t 2 ) s Vk,t r t 2 . 

(b) <f>(u ,t ) is a stationary random process with the function 

E{exp(j[<f>(u,t+x)-<j>(u,t)])} = *(t) 
known . 

(c) {A k (u,t)} and 4>{u,t) are independent random processes. 

Under these assumptions it is easily shown that 

" E{m.(u,t 1 )m k (u,t 2 )} 

-2 , Jw c^ d i" d k^ 

A^^-t^.+d^e c 1 K , i7k 

[K A ( t 1 - t 2) + A 2 ]$( t 1 -t 2 ), i-k 
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Defining Fourier transforms 


S m _ (f) - f R m ( T )e' j2 ' fT dT 

m i ra k J m iV 

—CO 


S^(f) '= j *(r)e 


~ j 2irf t. 


S A (f) = / K A (r)e 


-j2irfr 


It follows that 


2 J2*(f -f)(d--d k ) 

Irs (f)e c 1 , i7k 

$m i m k (f) ” 2 

1 k [A 6(f)+S A (f)]*S (f), i=k 


This can now be used to simplify the expression for the inter- 
ference power spectral density at the RF terminals of the j — 


subarray receiver, 


s j( f ) - 2N 0 |G RF (f)| 2 +[S A (f)*S (f)] £|G .( f ) 


-2 li ■ j2,ld k( f r- f ) 2 

+ A V f) Z V f)e 



This has considerably simplified the requirements for power spectral 
density information, and indicates the precise dependence of the 
interference power spectral density on the baseband equivalent RF 
system functions and on the pilot signals direction of arrival (w hich 
is hidden in {d^}). 


eJlinC^oi 


inborn. 


6.0 DESIGN PARAMETERS AND CONSTANTS 

The preceding sections provide a mathematical framework in which 
to study the interference rejection problem. We will now choose some 
design parameters and design constants for use in later computations: 

(a) Thermal Noise Density: 

N Q = kT = .69 x 10~ 20 

Here k is Boltzmann's Constant and T is 500°K. 

(b) Spacetenna Subarray Transmitter Noise Characteristics: 

1 

A = 1 v?P e "* 9 


2cx 

; S A (f) = 4Pe \s 0 (f) 


where 


P = Transmitter Power = 65 kw 
2 

= Transmitter Phase Error Variance 

3 • 

= .0305 rad 2 •v (10° rms) 

« = Variance to-squared-mean ratio of the angular modulation 

* 4 

S Q (f) = normalized phase error spectral density 

‘c , |f| < 1 KHz 

c(f/10 3 )" 6 , 1 KHz < |f| < 10 KHz 

U 10“ 6 (f/10 4 )~ 2 , 10 KHz < Jf I 

3 -1 

c = normalizing constant = (2.4 x 10 ) 

The development of the model for subarray transmitter noise 
is discussed in Appendix V. S ft (f) in above 1 W/Hz 
is shown in Figure 6.1. 
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(c) Spectral density of spurious phase modulation on the 
uplink signal : 

S (f) = (Dirac Delta Function) 

This is a reasonable approximation for interference 
computations. 


(d) Chip rate: 


T = 10 sec 


(e) Baseband-equivalent RF filter: 


- 


l+j2irfT 


The time constant t of this filter is normalized in terms 
of chip times to 

A T 

Ct *r 

c 

Note the 3 dB cutoff point for G^p(f) is tj— • Figure 6.2 i; 
a sketch of i G R p(f)| • The one-pole Butterworth filter 
characteristic has been chosen to conform with frequency 
response of half module 1 radiator supplied by Boeing. A 
comparison is given in Fig. 6.3. 

'(f) The noncoherent-interference-coupl ing equivalent-baseband 
system function: 

H i s kj( f) i 2 - K i ( constant ) 

k 

Here K-j is an array design parameter. 
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Figure 6.2. The Assumed Power Response of the RF Filter. 



AMPLITUDE 


Figure 6.3. Frequency Response of Half Module Radiator 
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(g) The coherent-interference-coupling equivalent-baseband 
system function: 


E G kj (f)e 


J2- IT d k (f c -f) 


^ (A Constant) 


Here K 2 is an array design parameter. 

(h) The PN apparent tracking offset e is normalized in terms 
of chip times. 


3 T. 


_e 

c 


(i) The first IF filter in equivalent baseband form is given by 


G i (f) " 1 ' + "ft tt ~ 

(jf) - exp(j(2n+l), T /12) 

This is a high-pass Butterworth filter with an added constant 
Y which can be used to specify the quality of the filter's null 
at bandcenter. Six poles were chosen since this matches the 
predicted theoretically optimal filter (see Section 4) in 
the critical range where S^(f) has a 60 dB per decade slope. 

The inversion frequency fg is a design parameter. A sketch of 
|G (f) | is shown in Figure 6.4. 


We ‘also normalize f Q to the chip rate by defining 


6 


f 0 T C 


(j) Bi-$ Manchester waveform characterization: 
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7.0 INTERFERENCE LEVELS AT THE OUTPUT OF THE SECOND IF FILTER 

Under the assumptions of the previous section, the interference 
spectral density in the worst case is given by 

Sj(f) - 2M 0 |6 RF (f)| 2 + K]S A (f) + K 2 s%(f) 

Using this as our model for $j(f), we can evaluate the effect of 
the notch filter parameters fg and y on the interference power Ip 
at the output of the first IF filter. 

h = / Sj(f) IG t ( f) | 2 df 

—CO 

The results of this computation are given in Table 1 7-1, along with 
a breakdown of the contributions to Ip of the three components of 
Sj{f), namely thermal noise, noncoherent downlink interference, and 
coherent downlink interference. The width fg of the notch filter 
is limited by the loss of uplink signal power which it causes. This 
will be evaluated in the next section. The minimum value of y is 
limited by hardware implementation problems. 

The effect of the PN code multiplier in~ the receiver is to spread 
the interference power spectrum over a wide band. The power spectral 
density of the input to the baseband equivalent second IF filter 
G^Cf) ,is given by 

3 I21N ( f ) - S c (f )*[ |G-, (f ) j 2 Sj(f)] 

R 

.where for an m-sequence, the code spectral density consists of lines 
1/MT C apart with: 

s c = ^ P n 6 D^ f 'MT"^ 

R n c 
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Table 7-1. Effect of Notch Filter Parameters On Interference Level Before 
PN Despread. 

(a) Effect of Cutoff Frequency 


f 0 ' 

X F 

I?"” 

6.94(1) 

3 

10' 3 

1.35(1) 

5 x 10 3 

7 . 1 0( — 3) 

10 4 

1 .90(-3) 

10 5 

1 . 30( -3 ) 

10 6 

1 .30(-3) 

2 x 10 6 

1 . 30( -3 ) 

3 x 10 6 

1.30(-3) 


COMPONENTS OF It 


HtRMAL 


7 . 83 ( - 1 3 ) 
7.83(-13) 
7. 82 (-13) 
7. 82 (-13) 
7.80(-l3) 
7.54(-13) 
7.26(-l 3) 
6 . 99 ( — 1 3 ) 


NONCOHERENT 


6.94(1) 
1.35(1) 
5.88(-3) 
6.78(-4) 
7.57(-5) 
7 - 57 ( —5 ) 
7 . 57 ( -5 ) 
7. 57 (-5) 


COHERENT 


1.22(-3) 

1.22(-3) 

1 . 22 (—3 ) 

1 . 22 ( -3) 
* 

1 „22(-3) 
1.22(-3) 
1.22(-3) 
1 .22(-3) 


(b) Effect of DC Attenuation with Narrowband Notch Filter 


COMPONENTS OF I, 


X F 


1.30{-1) 
1 .90(-3) 
6 . 1 7 ( -4) 
6.04(-4) 


' cJ^i i 
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(c) Effect of DC Attenuation with Wideband Notch Filter 


f 0 = 2 x 10 6 , B = o, a = , K 1 = K 2 = 10~ 2 , T c = 10" 7 



COMPONENTS OF I p 

J F 

THERMAL 

NONCOHERENT 

COHERENT 

1.30(-1) 

7.36(-13) 

7.57 ( —3) 

1 - 22 ( — 1 ) 

1 .30(-3) 

7 .26 (-13) 

7-57 (—5) 

1 ,22(-3) 

1.30(-5) 

7 . 26 ( -13) 

7.57(-7} 

1 .22(-5) 

1.30(-7) 

7 .25 (-13) 

7.57(-9) 

1 . 22( -7 ) 
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where 


1 1 

P E (0) 

7 

T c 


M ( 1+ l? 


p eW ^ 


, n = 0 


» n f 0, 


Hence the resultant spectral density at the input to the baseband- 
equivalent second IF filter contains lines separated by hertz 
along with a continuous spectral density. Vie assume that the 
second IF filter has a bandwidth narrow enough to reject all lines 
except the center line in ' For exam P le let’s assume that G 2 (f) 

is a zonal filter: 


¥ f > - 


{ 1 ’ 
( 0 , 


Ifl 


2MT, 


otherwise . 


This idealistic assumption is not critical since a "narrower" non- 
ideal filters occur in the carrier tracking loops which follow. Then 
we approximate the output baseband equivalent second IF filter by: 

S I20l!T^ = 2N eq + 2J( V f ) ’ < 2MT - 

c 

where 2N is the continuous spectral density at band center, and J 
is the d.c. line power, 

2 2 

2N eq ~ S P JVmtT^ [ 2N 0 ^RF^MTT^ + K 1 S aWH] 5 
TIC c c 


2d = y 2 P 0 K 2 S 2 

■Notice that if P £ (t) is a perfect Manchester pulse, i.e., e = 0, then 
J = O.and the line in this density disappears, implying that in this 
case the interference can be made proportionately smaller by decreasing 
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the loop bandwidth. 

The spectral density s j2out^ whiich We have J’ ust developed is 
the spectral density of the complex valued modulation on the IF signal 
at the output of Hg(f). The spectral density of the actual IF signal 
at the output of Hg (f ) is given by 

N 

IF y 


^+T[«Cf-f TP ) + <s(f+f IF )] } |f±f TF l < 1 


2HT 


S Sj20UT^ f ^ 


0, 


otherwise 


This expression will be used in later sections to determine tracking 
loop signal -to-noise ratios. The effect of the code length M on the 
noise process is shown in Table 7-2. 

8.0 UPLINK SIGNAL LEVELS AT THE OUTPUT OF THE SECOND IF FILTER 

The uplink pilot signal power P R at the terminals of the subarray 
antenna is given by 

_ p t G T A R 


4nR 


where 


Pilot transmitter power = 65 kw 
4ttAt 

n — -J- ~ Pilot Antenna Gain = 32,900 


= Pilot antenna aperture area = 78.5 nr (10m dish) 
= Transmit efficiency = .5 


P T 

g t = 

A 

A = Pilot transmitter wavelength = .12245m 
n 

f 

A 0 - Subarray pilot antenna perture = 3n^ 

K 

R = range = 38,000 km 


'cJClnC^c 
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The received pilot signal power then is 


P R = .354 yw. 

As discussed in Section 3.1, the pilot signal undergoes a 
fractional gain due to the notch in the first IF filter and to RF 
filtering, which is given by 


L 


F 



|P E ( f )| 2 G RF (f)G l (f)e " i2Trfedf 


where e is the PN apparent tracking offset. FigureSrl indicates the 
variation in Lp as the normalized PN apparent tracking offset 3 varies 
between _+! chip time. Notice that in the most optimistic case for the 
chosen parameters, Lp is less than 0.7 and that tracking error from this 
maximum must be on the order of 0.1 T to keep L- above 0.6. Our nominal 

t C r 

choice of g = 0 for later power computations puts Lp at about 0.55. 

The fractional power loss is also function of the RF terminal response 
G R p(f) and the notch filter G-j(f). The results are given in Table 8-1. 
For these tables, we have also indicated the loss if phase compensation 
for the filters (G R p(f) and G-j (f ) ) are used. 

We assume that the carrier tracking loop is operating properly and 
that a nominal unit gain is seen by the pilot signal's center-line when 
passing through the second IF filter. 

9.0 COSTAS LOOP PHASE ERROR 

At the input to the Costas loop, the signal consists of three terms: 

a term due to the uplink pilot data modulation with power P R *Lp, a line 

due to the interference with power J and an equivalent "white" noise term 

with one-sided spectral density N . From [7 ], the phase error of a 

* eq 

Costas loop is given by 





Table 8-1. 


Effect of G p(f) and Notch Filter Bandwidth on 
Fraction PoWer Loss Lp. 


(a) Effect of Notch Filter 


6 

L F (peak) 

* 1 * 

L p 

10" 3 

.63 

.63 

10" 2 

.63 

.63 

10" 1 

.63 

.63 

.2 

.51 

.61 

.4 

.35 

.51 

.6 

.29 

.35 

.8 

.15 

.18 




a = T 


(b) 


Fffect of RF Terminal Response G 


RF 


(f) 


(2ira) 1 

1 * 
l f 

.2 

.06 

.5 

.25 

1 

.49 

2 

.72 

5 

.88 

10 

.93 

20 

.96 

50 

.97 


T c = 10- 7 

Y = 10“ 3 

6 = TO" 1 


*With ideal phase compensation. 
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-i 

L V l f 


where is the one-sided loop bandwidth and ^ is the squaring loss 

which is a function of the IF filter and the arm filters. A typical 

squaring loss is 2 dB. The interference term will also introduce a 
' \ 
phase error which depends on J and the random phase between the 

interference and the carrier. We shall adopt a worst case analysis, 

i.e., the carrier and the interference are at quadrature with each 

other. The phase error is now modified to 

-1 % B L , J 

L p r' l f p r’ l f 


If no data modulation is employed for the uplink pilot, the carrier 
recovery can be accomplished with a CW loop. In that case, the squaring 
loss penalty is removed, it also reduces hardware problems such as balancing 
the I and Q arm filters when a Costas loop is used. 

10.0 SPREAD SPECTRUM CODE SELECTION 

The system design, as analyzed earlier, depends on the spread 
spectrum sequence (SSS) {a^ } in several ways: 

(a) The d.c. line energy in the code spectral density, namely 
Pq, is the same for any balanced sequence of odd period. 

Hence we have assumed that the SSS has a number of ones 
and a number of zeros in one period which differ by 1. 

(b) The SSS spectral line spacing is (MT C ) -hertz and it is 
assumed that all but the d.c. line would be rejected by 
the second IF filter. Hence the code period M must satisfy 
the relation 
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M « B 9 T 

2 C 

where is the bandwidth of the second IF filter, or the 
filter must have carefully constructed nulls at multiples 
of 1/MT C . 

(c) The equivalent noise spectral density parameter N is 

eq 

a function of the SSS spectral line levels P n for all n. 

(d) To provide a good characteristic for the code tracking 
loop, the SSS must have a spike-like periodic autocorrelation 
function. 

Now lets investigate the options open to us under these constraints. 

With a chip time of lO - ' 7 seconds and a second IF filter bandwidth 
on the order of 100 hertz, the period of the sequence must be less than 
10 to satisfy the line rejection constraint (b). Certainly maximum- 
length linear-feedback shift register sequences; i.e., m-sequences , 
satisfy the remaining constraints. They are balanced and have spike- 
like autocorrelation functions. The numerical data for N given in 
Section 7 is based on the values of P n for m-sequences. 

The following gives the number of distinct m-sequences and their 
period M, as a function of shift register length L, along with a 
typical shift register generator. 


Register Length 

Period 

Number 
of Sequences 

Typical 

Generator 

9 

511 

56 

9,4 

10 

1023 

99 

10,3 

11 

2047 

186 

11,1 

12 

4095 

335 

12,6,4,1 

13 

8191 

630 

13,4,3,1 

14 

16383 

1161 

14,12,2,1 
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To construct the m-sequence from the generator information in the 
above table, simply take the modulo 2 sum of the outputs of the stages 
indicated in the table and feed the result back into the shift register. 
An example is given below for the length 12. 

STAGE NUMBER 



Load the register with any sequence of bits, not all zero, and as it 
runs it will produce some phase shift of the m-sequence associated 

N 

with the particular register. Changing to a different seqeunce involves 
changing the tap connections. Appropriate connections for other 
generators are determined using a table of primitive polynomials over 
GF(2), e.g., the table in [3], 

If code-division multiple access schemes are required to provide 
satisfactory isolation between the pilot signals from different rectennas, 
then there are several alternatives, all based on shift register 
generated sequences. They are (a) Gold Codes, (b) Kasami Sequences, 
and (c) Bent-Function Sequences. All are easily implemented. Properties 
of these sequence sets are tabulated below for designs of period 2 n -l . 

The Gold Codes are suboptimal from two viewpoints: They are not all 
balanced, and their crosscorrelation is approximately 72 too large. 

Both Kasami sequences and bent sequences offer an alternative to Gold 
sequences with better correlation properties. The Kasami sequences 
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Property 

Gold [4] 

Kasami [5] 

Bent [6] 

All Balanced 

No 

No 

Yes 

Number of Sequences 

2 n +l 

2 n/2 

2 n/2 

Allowable n 

n not divisible 
by 4 

n even 

n a multiple 
of 4 

Maximum Cross - 
Correlation 

2 (n/2)+l +ljn even 
2 {n+l)/2 +1 ^ n odd 

2 n/2 +l 

2 n / 2 + i 

Linear Recursion 
Length 

2n+l 

3„ 

2" +1 

Longer 


are available in a greater variety of lengths (2 n -l for all even n) than 
bent sequences, but the bent sequences are guaranteed balanced, and in 
addition have two nice properties when an intelligent jammer is trying to 
cause trouble. First, the linear span of a bent sequence is quite long 
for the shift register length employed. Secondly, it is quite easy in 
hardware to reinitialize a given bent sequence to a new random phase. 

In the code division multiple access situation, we would recommend 
the use of the set of 64 bent function sequences of period 4095. The 
design details are available in [6] along with a sample design for n=12. 
For a further discussion, see Appendix 6. 
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11 .0 PERFORMANCE EVALUATION VIA SOLARSIM 

In the precding sections, analytical models have been developed 
in terms of the pertinent design parameters of the receiver portion 
of the SPS transponder. These important parameters include; 

©Bi-4 Waveform Asymmetry 

©PN Chip Rate 

©PN Code Period 

©Uplink EIRP 

©Receiver G/T 

©Diplexer Isolation 

•RF Path Cutoff Frequency 

eNotcn Filter Stop Band and Attenuation 

©Costas Loop Bandwidth 

They serve to characterize the uplink waveform, the Jn'nk budget, the 
receiver RF front end, interference nulling IF filter and the-'phase 
tracking system. In order to optimally select these parameters, the 
SOLARSIM has been updated to include a computer program to evaluate 
the performance of the transponder as a function of these inputs. The 
performance measure selected is the phase error of the Costas tracking 
loop. A description of the computer package is included in Appendix 6. 

To demonstrate its value as an evaluation tool, we have used the 
SOLARSIM package to study the effect of (1) cutoff frequency of RF 
system function, (2) notch filter cutoff frequency, (3) notch filter 
attenuation and (4) PN code period on the Costas loop tracking phase 
jitter with the PN chip rate as a variable. For these examples, the 
nominal conditions are chosen to represent ,a practical transponder 
design and they are listed in Table 11-1. The simulation results are 
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plotted in Figs. 11.1 - 11.5. 

Figure 11.1 indicates that widening the RF filter bandwidth relative 
to the chip rate improves the system performance. However, the improve- 
ment margin becomes smaller as f R p/R c (f^p is the 3 dB cutoff frequency) 

* 

increases, indicating there is a region of diminishing return beyond 

f^p/R c = 2. This is to be expected since most of the pilot signal energy 

is concentrated for f < 2R c> In the present system, the RF bandwidth 

is pretty much set by the antenna frequency response at about 20 MHz. 

A good choice for the chip rate is therefore around 10 M chips/sec. 

In the next figure, we see that the phase jitter reduces inversely 

with the notch filter bandwidth for a fixed chip rate. This phenomenon 

continues up to f^/R c ^0.5 (f^ is the 3 dB cutoff frequency for the 

notch filter) and the trend reverses itself as indicated in Fig. 11.3. 

In the first region, the notch filter suppresses the interference v/ith 

a minimum distortion on the signal. After a critical point is reached 

the notch filter suppresses both the interference and the signal, and 

penalizes the latter more severely. The effect of the band center 

attenuation of the notch filter is. shown in Figure 11.4. As to be 

expected, higher attenuation improves phase jitter performance. 

For R c = 10 M chips/s, the phase noise contribution (suppressed) 

becomes negligible and the performance does not improve appreciably 
2 

with increasing y . In general, since the cost of the notch filter goes 
up with the attenuation, it seems logical to choose a value that is 
cost effective yet meeting the performance requirement. 

From Figure 11.5, we can observe that the phase jitter performance 
improves as the PN code period M is increased. As we have discussed 
earlier H is upper bounded by R c /Bp. A reasonable choice is M % 10 4 . 
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in L/om 



CHIP RATE (HZ) 


79 0269 
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12.0 REFERENCE SYSTEM SPS POWER TRANSPONDER 

So far, we have considered the pilot receiver portion of the 
power transponder in conjunction with the selection of the uplink 
pilot signal parameters and the design of the transponder RF front 
end. In the following sections, we shall evaluate the overall system 
performance of the transponder. In particular, we shall proceed to 
characterize the transponder's ability to perform the phase conjugation 
in terms of the key technical issues specific to the SPS environment. 
Figure 12.1 represents the overall functional diagram of the" SPS power 
transponder. This includes the pilot signal receiver, phase conjuga- 
tion electronics and the high power amplifier phase control subsystem. 

In the mechanization of the SPS power transponders, two 
receiver "types" will be required; however, most of the hardware 
will be common between two receivers. One receiver, the Pilot Spread 
Spectrum Receiver , is located at the center of the spacetenna or the 
reference subarray. It serves two major functions; (1) Acquires 
the SS code, the carrier and demodulates the command signal, (2) 
provides the main input Signal to- the Reference Phase Distribution 
System, see Figure 12.1. 

The second receiver "type" will be located in the Beam Forming 
and Microwave Power Generating System, see Figure 12.2. Its main 
purpose is to phase conjugate the received pilot signal and trans- 
pond power via the j-th spacetenna element, j = 1 ,2, . . . ,101 ,552. We 
now discuss the functional diagram indicating the mechanization of 
the SS power transponder and discuss its operation. 

From Fig. 12.1 we note that each SS receiver -element must be 
capable of despreading and demodulating the received pilot signal 
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Figure 12.1. Central SPS Power Transponder Located at Spacetenna Center. 
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based upon the gain provided by a single element of the SPS space- 
tenna. The receiver consists of several major subsystems. These 
include, (see Figs. 5.1 and 5.2): (1) the RF Front End centered at 
2450 MHz, (2) The SS Code (PN) Sync Subsystem (PNSS), (3) The AGC 
Subsystems, (4) The Carrier Sync Subsystem (CSS), (5) The Carrier Lock 
Detection Subsystem, (6) Carrier Sync Acquisition Subsystem, (7) Symbol 
Synchronization Subsystem, (8) the Viterbi Decoder (if the uplink 
employs convolutional coding), and (9) The Ambiguity Resolving Subsystem. 
Note that (8) and ‘(9) are not required if the pilot signal is not data 
modulated. 

The carrier, with nominal frequency of 2450 MHz, is first processed 
via an RF filter. The bandwidth of this filter must be sufficiently 
wide to pass the PN chips and the roll-off must be sufficient to meet 
the desired rejection requirements. 

The signal level into the first IF mixer is held constant by the 

noncoherent AGC. The first L0 is selected to run at 1960 MHz; 

/ 

therefore, the first IF frequency, at .zero Doppler, is 490 MHz. The 
output of the first IF mixer is further filtered by the first IF filter 
whose bandwidth is sufficiently wide to pass the PN chips. The cascaded 
frequency response of the RF filter and the first IF filter are collect- 
ively designed to meet the desired front end rejection requirements. 

The IF filter output serves as the input to the PN synchronization 
system. The PN sysnchronization system ( J PNSS) of Figure 12.1 incor- 
porates a noncoherent x-dither PN acquisition and tracking design. 

The arriving NRZ/BPSK/BI-<j>-DS/CDMA signal is despread prior to 
filtering by the second IF filter. This gives rise to an ordinary 
BPSK signal when the x-dither loop is locked. The output from the 
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second IF filter goes to the output of the carrier recovery circuit. 

The bandwidth of the second IF filter is chosen to be wide with 
respect to the data rate. 

A Costas (I/Q) loop configuration is chosen for carrier 
acquisition, tracking and data demodulation. This configuration 
was chosen because it was determined to be optimum when all 
considerations, including the ability to square perfectly over 
temperature and signal level, are traded against lock detection arid 
synchronization monitoring. 

Noncoherent AGC is derived from the sum "I" and "Q" Channels 
appearing in. the arms of the Costas loop. In addition, lock detection 
for the carrier circuit is accomplished by using the difference between 
the "I" and "Q" channels of the Costas arms. The noncoherent AGC£ 
controls the receiver gain (prior to phase detection} with the signal 
plus noise level appearing in the outputs of the Costas loop arm. This 
feature is used to control the loop bandwidth and damping factor during 
acquisition and tracking. 

The loop filter receives the signals from the phase detector 
(third multiplier) and supplies an error signal to the VCO which 
control the local frequency. The loop filter sets the tracking 
loop bandwidth and damping factor. An AGC voltage, proportional 
to the incoming signal plus noise power, is low pass filtered and 
amplified to drive variable gain elements in the IF amplifiers. 

(These amplifiers are included here in the first and second IF filters 
for simplicity.) Additional integrate and dump circuits and threshold 
detector circuits control the sweep and the lock indicator needed for 
loop supervisory control. It is to be noted that the carrier sweep 
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is not activated until PNSS is synchronized. 

The carrier lock detector circuit monitors the integrate and 
dump voltage formed at the discrete points in time by differencing 
the squares of the inphase and quadrature arm outputs. .The output 
of the integrate and dump circuit is compared to a fixed threshold 
level to detect lock. 

If the pilot signal is data modulated, the data can be extracted 
from the output of the Q-channel as shown in Fig, 1 2 . 2. The 180 degree 
phase ambiguity the Costas loop introduced can be resolved by oeriodically 
inserting a (predetermined) identification sequence in the transmitted 
data stream; the recovered data can be checked against this known 
pattern to correct for the phase ambiguity. When the pilot signal is 
not data modulated, the Costas loop can be replaced by a CW loop. 

This avoids the need for provisions in the carrier tracking subsystem 
to resolve phase ambiguity. 

The 490 MHz reference appearing at the pilot receiver VCO output, 
see Figure 12.1, serves as the only input to the Phase Control System. 

This phase characteristic is then distributed over the aperture of the 
spacetenna. The output of the Phase Control System consists of 101,552 
980 MHz constant phase reference signals which are used to conjugate 
the incoming signal and to stabilize the high power amplifier (HPA) 
outputs. This is achieved by placement of an automatic phase control 
around each HPA, see Figure 12.2. 

13.0 SIGNAL- AND NOISE CHARACTERISTICS 

For the purpose of the SPS transponder analysis, the signal and 
noise spectrum into the transponder can be represented as shown in 
Fig. 13.1. The spectral shape of the input Bi-$-PN spread pilot signal 
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Figure 13.1. Signal and Noise Spectrum into SPS Transponder. 
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4 2 

is of the form sin (x)/x with the nulls occurring at frequency 
multiples of 2R (R = PN chip rate). The receiver front end 
contributes an equivalent thermal noise level equal to W/Hz. 

The phase noise interferences from the power beams (including its 
own) can be modeled by the coherent and noncoherent terms and indicated. 
The important RF filter and notch filter characteristics are also shown 
in the same figure. The RF filter is used to model the bandwidth of 
the receiving antenna (waveguide arrays). The notch filter is intro- 
duced to suppress the phase noise interferences. It is clear that 
since the portion of the signal power close to the carrier is small, 
the signal .distortion introduced by the notch filter is negligible 
while the interferences can be significantly suppressed. The optimal 
selection of the filter parameters have been considered in the preced- 
ing sections. 

The klystron amplifier in the power transponder generates a phase 
noise profile that can add significantly to the downlink carrier phase 
error if it is unattended for. Figure 13.2 shows the phase noise 
sideband power spectral density measured on a Varian X-l 3 klystron 
amplifier,* which we have adopted as our klystron model. Upon 
reflecting to 2.45 GHz, the phase noise component that is outside a 
bandwidth is plotted in Fig. 13.3. Selected values are tabulated 
in Table 13-1. The components of the klystron phase noise around 
the carrier frequency can be tracked with the PA phase control loop 

a 


provided for in the reference system. With the loop around the klystron, 
only phase noise components which have Fourier frequencies greater than 


*The complimentary measurement is provided to Lin Com by Dr. Algie Lance 
of TRW, Redondo Beach, CA. 
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will be transmitted. Components below this frequency will be 
suppressed by a factor l-H(s) where H { s } is the closed-loop transfer 
function of the PA control loop. Assuming an ideal loop transfer 
function, namely H(f) = 1 for f < and H(f) = 0 for f >_ B^, we can 
predict the amount of phase noise leakage through the transmitter. 

For example, from Table 13-1 a 400 Hz loop will introduce a rms phase 
error of 4 degrees. When the phase noise effect is the only consider- 
ation, the obvious solution is to widen the loop bandwidth as far as 
possible. However, this also opens up the bandwidth for other noise 
contributions through the transponder. Clearly, a tradeoff has to be made 
to yield an optimum balance. A more detailed modeling of the overall 
system which enables such a tradeoff will be addressed in Section 14.0. 

Other noise sources that affect the performance of the transponder 
includes the VCXOs, mixers', multipliers, up converters and down- 
converters. Their relevance will also be discussed in Section 14.0. 


14.0 SPS TRANSPONDER TRACKING LOOP SUBSYSTEMS 

At the proposed synchronous orbit, the channel Doppler profile 

is characterized by a maximum Doppler of 25 Hz and a maximum Doppler 
-3 

rate of 2 x 10 Hz/sec, at the carrier frequency of 2,450 MHz. For 
this reason, loop acquisition does not pose any significant problem. 

The loop bandwidths can be made to be relatively small (as low as two Hz) 
-and still guarantees satisfactory response time. 

14.1 PN Tracking Loop 

The main purpose of the PN tracking loop is to despread the 
incoming Bi-<f>-PN signal. It is not designed for .Doppler or range 
measurement. As such, the tracking error requirement is not very 
critical as long as it is limited to a few percent of a chip time. 
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PN acquisition will take time from a "cold start," i.e., when 
the local PN clock is not running. This can occur, for example, when 
the system is powered on initially. In that case, the local clock has 
to step over all possible code states during acquisition. However, 
once the local clock starts running, it should be almost synchronous 
with the ground code due to the almost nonexistent channel Doppler. 

If the local clock is left running when the system is brought down 
for whatever reasons, the latter can be brought back up using the 
local clock. No additional acquisition algorithm is required. 

Since all power transponders are experiencing the same Doppler 
for all practical purposes, a way to cut down the acquisition time 
and individual acquisition hardware requirement is to include a separate 
telemetry receiver that tracks the uplink pilot. This pilot signal is 
constantly tracked by the telemetry receiver and the state of the 
local PN clock can be transferred to the individual transponders to 
start the PN loops. In that case, no acquisition aid on the trans- 
ponders are required. However, data links between the telemetry 
receiver and the transponders must be established. 

In the reference system, a -r-dither loop is used to avoid the 
gain imbalance problem commonly found in the standard "delay-lock loop 
implementation. Fig. 14.1 gives an alternate implementation [8] of 
the delay-lock loop which is immune to gain imbalances and is equivalent 
to the classical delay-lock loop. The delay-lock loop can potentially 
outperform the T-dither loop by approximately 1 dB. 

14.2 Phase Reference Tracking Loop Model 

Fig. 14.2 shows the equivalent mathematical model of the phase 
reference tracking loop located at the center of the spacetenna. 
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The phase reference tracking loop is mechanized as a long loop. The 
statistics of the loop phase error ^ and the phase of the reference 
e R (both measured at 490 MHz) are of importance to assess performance. 
In the tracking mode, they are related to the input noise processes via 


1-H-, (s) H ] (s) / 5n.j \ 

= ThFTjTsT ( 0 o~ 5, h ^ ■ T+4FTjTs7 \7TJ 

l-H-j (s) H-j (s) / n] \ 

s r ~ mFfjTsT + + Q oJ 


(14-D 


04-2) 


where 8 q is the transmit phase noise including the ionospheric disturb- 
ances, n^ is the equivalent thermal noise seen by the loop, P is the 
transmitted power seen by the loop, models the equivalent VCO/mixer 
phase noise and H-j (s ) represents the closed loop transfer function of 
the Costas loop. The squaring loop mechanization also introduces a 
squaring loss wnich is a function of the arm filters and the IF filter 
bandwidth preceding the I-Q phase detectors. This factor must be 
accounted for, presumably in computing n-j before using (14-1) and (14-2) 
for performance evaluations. 

14.3 Power Transponder Phasing System Model 

The equivalent mathematical model of the individual transponders 
is given in Fig. 14.3. The relevant tracking equations for the pilot 
■phase recovery loop are 


«? 2 = [l-H 2 (s)][6(.e R +<j) D )-e 1 ] - [1-H 2 (s)> 2 - H 2 (s) (14-3) 


and 


e 2 - [1-H 2 (s)]4i 2 + h 2 (s) 


7F 


+ 6(e R +* D ) - e ] j 


(14-4) 


where «$ 2 is the tracking error, the transmit pilot phase noise, ^|> 2 
models the equivalent VCO/mixer phase noise, n 2 is the equivalent thermal 
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noise, P the transmitted pilot power, (s ) is the closed loop transfer 
function for the Costas loop and <j>p models the phase noise and 
differential delay introduced by the phase distribution system. The 
relevant tracking equations for the PA phase control loop are 


% = D-H 3 (s)]9 2 - [1-H 3 (s) 3(+ a +* 3 ) (14-5) 


e out = Cl -H 3 (s )] (<j > A +^ 3 ) + H 3 (s) 0 2 + 4(e k +<j> D ) (14-6) 

where «P is the loop phase error, H 3 (s) is the closed loop transfer 
function, is the (PA) klystron phase noise, and ^ models the VCO/ 
mixer phase noise. The noise process e oLjt is the single most important 
quantity as it models the phase error process at the output of the 
transponder and it directly affects the efficiency performance of the 
SPS. Note that except for the klystron phase noise which is measured 
at 2450 MHz, all other phase noise processes are measured at 490 MHz. 

14.4 Overall Transponder Equivalent System Model for Analysis 

The analytical model for the performance of the overall transponder 
is given in Fig. 14.4. Notice that the inputs to the transponders are 
0 q, 0-j , n-j apd n 2 ; 6 q and 0-j represent the phase noise on the pilot 
signal including the effects of the ionospheric disturbances and transmit 
frequency instability and n-j and n 2 represents the effects of the receiver 
thermal noise, power beam interferences ('caused by mutual coupling) and data 
distortions. The phase disturbance at the output of the transponder is 
the term of interest as it directly affects power transfer efficiency 
of the SPS. We have studied this effect extensively using the SOLARS IM 
program. All other noise components are generated within the power transponder 
They are the filtered version of various hardware induced noises. 
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Figure 14.4. 


Overall Transponder Equivalent System Model 


for Analysis. 
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In Fig. 14.4, filters of the form H(s) are low pass filters and 
filters of the form l-H(s) are high pass filters. Note that oscillator 
phase noises (including the klystrons) are all filtered by a high pass 
filter. Since oscillator phae noise profiles are typically of the form 
indicated in Fig. 14.5 with most of the phase noise power concentrated 
around the carrier, one should widen H(s) as far as possible. For 
example with our klystron model, the bandwidth of H^s) should be 
kept to approximately 1 KHz to limit the phase noise introduced jitter 
to less than 1 degree. However, widening H(s) also introduces other 
problems. As another example, opening up H 2 (s) and H^(s) allows 
more phase jitter due to n 2 and e-j. Obviously, in order to attain the 
optimum choice of loop bandwidths, we need to specify power spectral 
densities of all the noise processes, and then perform tradeoffs based 
on the system model. In view of the fact that many parameters are 
involved, a simulation program is needed to exercise these tradeoffs. 

A careful examination of Fig. 14.4 also reveals the following 
important fact: some of the noise processes are irreducible. For 

example, one can reduce the effects of n 2 //P*by changing the system 
designs such as reducing the bandwidths of H 2 (s) and H 3 (s), increasing 
the uplink power", . increasing the processing gain of the PN signal 
or employing diplexers with improved isolation characteristic, etc. 
However, there is nothing one can do about the random phase disturbance 
introduced by the phase distribution tree because a major portion 
of appears at the transponder output unmodified. From Fig. 14.4, 
one can see that the irreducible error terms are <j>g and the ionospheric 
differential phase perturbations e Q and . -(The processes e Q and e-j 
can have a significant portion of power within typical values of loop 
bandwi dths . ) ^ _ 

csLinL^om 
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15.0 RECOMMENDATION FOR OVERALL TRANSPONDER DESIGN 

Although a thorough simulation study is required to quantify 
various tradeoffs, one can still make some, general recommendations 
on the transponder design. From the power transfer efficiency study, 
we know that the transponder phasing system have to hold its rms phase 
error to under 10 degrees. In that case, it is reasonable to limit 
the contributions from each individual error source to less than 1 
degree rms. Under this guideline, we can make the following 
recommendations: 

(1) Klystron phase control loop bandwidth >_ 10 KHz. 

(2) Code loop jitter can be made negligible and is not critical 
to performance. 

(3) . Costas loop bandwidth fJ- 10 Hz. 

When these conditions are met, it is reasonable to expect that' the rms 
phase error introduced by the transponder, interference and thermal 
can be limited to less than one degree. However,, this does not take 
into account the effects of the reference phase distribution error and 
the hardware induced differential delay variations that differs from 
transponder to transponder. Unfortunately, these effects belong tothe class 
of irreducible errors that cannot be controlled by the transponder 
•design . 
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APPENDIX 1 

BASEBAND-EQUIVALENT FILTERING 


Let 


' ( Joj t| 

x(t) = Rejm(t)e c j 

JU3 t , -ju_t 

ra(t)e c + j m*(t)e 


- 1 
2 


Taking Fourier transforms, 

X(f) » ^H(f-f c ) +j-M*(-f-f c ) 

Suppose we pass this through a filter with real impulse response 
H(f) = H*(-f). Then the output y(t) has transform 

Y(f) = H{f)[lM(f-f c ) +lM(-f-f c )] 

Now if |M(f ) j = 0 for Jf| >_ f , then 

Y(f) = \ U(f)H(f)M(f-f c )+ |u(-f)H*i-f)M*(-f-f c ) 

Hence since Y(f) has conjugate symmetry 

y(t) = Re j J H{f)M(f~f c )e jblt df J 
0 

- Reje 3 " c j U(f+f c )H(f+f c )M(f)e J "“ t df | 


y(t) can now be modelled .as follows: 







APPENDIX 2 

BASEBAND-EQUIVALENT MULTIPLICATION 


Ocii-« t 

x(t) = Re{a(t)e } 

jcop't 

y(t) = Re{b(t}e } 


and suppose that 


z(t) = x(t)y(t) . 


■1 j +(JJo j (<dl -(Op )t 

z(t) = j Re{a(t)b(t)e + a(t)b*(t)e } 

Now suppose that a(t) and b{t) are bandlimited so that 

A(f) = 0, |f| > B a 

B(f ) = 0, |f| > B b 

Then the two components of z(t) can be recovered by filtering 

j 1 ^Re{a(t)b(t}e } 


o4?G{3(t ) b*(t)e 


f l > 

f 2 > ®A + B B » 0 

J 

1 , 

I|f|-(f 1+ f 2 )| < B a + b. 

0 , 

ll f l-<V f 2 )l < B A + B l 

1 , 

|| f |-(f 1 - f 2 i| < B a + B. 

0" > 

I if l-Cf,+f 2 ) 1 < b a .+ b. 




APPENDIX 3 

THE PRODUCT FILTERING APPROXIMATION 
Suppose that 

z(t) « x(t)y(t) 

and that z(t) is the input to a linear filter whose output is 
w(t), 

03 

w(t) = J x(a)y(a)h(t-a)da 

.CO 

Now suppose that x(t) is reasonably constant over the periods 
corresponding to the duration of the impulse response h(t). 
is then possible to approximate w(t) by 

CO 

w(t) £ X(t) f y(a)h{t-a)da 


This assumes h(t) is nonzero primarily near t=0. 
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APPENDIX 4 

SS CODE SPECTRAL DENSITY 


The SS code signal c R (u,t) is periodic with period MT C where 
M is the period of the sequence {a.. }. The correlation function can 
be found simply by time averaging the correlation product over one 


period 


V Tl ° MT C 


j_r 

W-t 


MT -T / 2 
c c ' 


c R (u,t+x)c R (u,t)dt 


with T (u) set to zero 


MT 


i M-i 

r e E E Vj J 

i=0 j -T /2 


MT c -T c /2 


P(t-iT c )P(t-jT c -x)dt 


Setting x - mT c +Xg, 0 £ Xg < T c , we can use the fact that P(t)=0 


for all |t| < T /2 to simplify this to 

V 


R (t) - R>)R p (t 0 ) + R a (ra+l)R p ( T 0 -T c ) 


R (in) = jt a. a. 

a v ‘ M Z-j i+m 

1=0 


R P^0> “ T 


P(t)P(t-Tg)dt 


■V 2 


It is apparent that R_ (x) is periodic with period MT and hence 

c R c 

has a line power spectral density 


s c ' f > ■ £ V (f - mtt) 


csCinC^ot 


where s{f) is the Dirac delta function and ? n is given by 




R (T)e- j2 * nT/MT Cd T 
Cn 



where p(f) is the Fourier transform of P(t) and S (n) is the discrete 

a 

Fourier transform of R (m): 

a 

S a (nj - 1 g R a (")e' ;i2 " n ' n/M 

m=0 


For the standard Manchester waveform 


p(f) = f P(t)e~ j27rft dt 


2 

-j-rrf 


sin 2 (7rfT c /2) 


whence 


4 sin^(7rfT_/2) 




APPENDIX 5 

THE SPECTRAL DENSITY OF A k (.u,t) 


The quantity A^Cu.t) is a complex-valued -multipl icative noise 
factor appearing in the downlink transmitter modulation for subarray k. 

It is assumed independent of A-{u,t), jYk. We shall suppress the 

J 

subscript k and expand A k (u,t) in an envelope and phase representation:, 

A(u,t) = lA(u J t)[e j0(u ’ t) 

Lets further assume that e(u,t) is Gaussian and independent of |A(u,t)| 
so that the autocorrelation function of A(u,t) is given by 

R A (t) - R |A| (T)R ph (T) 

where 

Rp h (t) = Ef exp( j[e ( u ,t+x )-e (u ,t)'] ) > 

Now e(u,t+ T )-e(u,t) is Gaussian with mean zero and variance 2af(l-p(x)). 

0 

2 

Here « is the variance of e(u,t). Using properties of the characteristic 
0 

function of a Gaussian random variable, it follows that 

R ph^ = ex P{“ a g 0 - P ( T ) ) } » 

where p( x ) is the normalized autocorrelation function of e(u,t). 

V/e assume that the specification on the downlink transmitter^ 
rms phase error is 10°, and that it is due mainly to transmitter 
generated oscillator noise. Then 


2 

a 8 " 


Notice further that 


/2tt\2 

W 


= .0305 rad 


c*Cin(^om 
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E{e de(u,t} ) = e a e /2 = 9Q 


Hence the covariance function of is given by 


K ph (x) = R ph (x) ~ e 6 

^0. V>( T ) 

= e 0 [e 9 -1] 


For our assumed value of phase variance, this is accurately approximated 


K ph (T) ^ 6 9 V (t) - 

The random variable exp{je(u,t)> has a variance-to-squared mean 


ratio a given by 


a - e 0 K . (0) = erf 


Assuming that |A(u,t}| has the same ratio a and the same normalized 
spectral shape as exp{je(u,t)} , it follows that' 

R a Ct) = R| A |(0)R^(t). 

Note that R ph (0) is unity and hence we define 

R a ( 0) = 2E{s 2 (t) } - 2P 


where P is the average transmitted power of a spacetenna subarray 


(P is nominally 65 kw), and 


2 


R a (t) = 2P[K .(t) + e °) 


~ 2Pe [ap(x) +1]* 


The mean A and covariance function K a (t) required in section 5 are 


-8 
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2 

- 


A = /2P e 0 


K a (t) % 4Pe- 0 ap(x) 

The only parameter which has yet to be specified in this model is the 
normalized correlation function of the phase noise e(u,t). 

Our computations actually use the Fourier transform of K^(x), 


namely S^(f). Hence 


2cj 

SAf) % 4Pe 00 «S n (f) 


where S Q (f) is the Fourier transform of p(x). Our model for this 


normalized density is 


s 0 (f) - 


C(f/10 3 )" 6 , 


[ f 1 < 1 KHz 
1 KHz < |f| <10 KHz 


V 10“ b C(f/10 4 )"\ 10 KHz < ff i 


where the normalizing factor C turns out to be 

C = (2.4 x 10 3 )' 1 . 


Jhis model is derived from the shape of the spectral density of a 
Varian X-13 kylstron tube indicated in Fig. A5-1. It is more pessimistic 
than the klystron data since more power is placed outside 1 KHz. 
(Presumably, the noise power inside 1 KHz will be notched out and 
does not affect performance.) 
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APPENDIX 6 
OPTIMAL CDm CODES 


This appendix surveys several known sequence sets which may be used 
as spread spectrum codes in code-division multiple-access communication 
systems. All of the designs basically achieve Welch's lower bound on 
the maximum value of periodic cross -correlation between signals and 
are optimal in this sense. 

1 . Introduction 

In 1974 L. R. Welch [1] published a bound on inner products which 
could be specialized to the case of periodic correlation of spread- 
spectrum code-division multiple-access (CDMA) signal sets. Specifically, 

j 

consider a set of M sequences a^., i = 1,...,M, of period L, 

a t = a t+L Vijt * (A6-1) 

The periodic cross-correlation between sequence i and sequence* j at 
shift x is defined as 


V t > 


& 


L-l 


E 

t=0 




(-A6-2) 


((*)* denotes conjugation), the maximum autocorrelation of the set is 


C, - max max |C. . (x)| 
‘ i 0 <t<L 11 


(A6-3) 


and the maximum crosscorrelation of the set is 


Co - max max |C.,.(x)|, (A6-4) 

c i&j 0 <t<L 

Under the assumption that the sequences all have the same "energy" 
per period, i.e.. 
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c it (°) = C..(0} VI J , (_A6— 5 ) 

Welch demonstrates the normalized correlation bound 


r ' a max f c r c 2^ /FTP 

max ~ C-j -j (0) - *} ML-1 


(A6-6) 


This bound has become the standard against which a possible CDMA 

signal set design is compared, despite the fact that C often 

max 

is not specifically a parameter in the communication system design [2], 

The purpose of this Appendix is to review several of the CDMA 
sequence designs which are known to nearly achieve the Welch bound. 

2. Optimal -CDMA Code Designs 

A. Kasami Sequences [3,4] were designed originally as a linear cyclic 
error-correcting codes. The underlying arithmetic in Kasami’s design is >i 
performed in the finite field 6F(2 n ), n even, with M-j(z) representing 
the minimum polynomial over GF(2) of a primitive element a of GF(2 n ), 
and r%'(z) representing the minimum polynomial over 6F(2) of c? , where 
s = 2 n ^+l . Thus a has order 2 n//2 -l and is a primitive element of 
GF(2 7 ). Hence M-j(z) and I^(z) can be viewed as the characteristic 
polynomials of binary (0,1) linear feedback shift registers which 
generate m-sequences, b| and b 2 of lengths 2 n -l and 2 n ^-l respectively. 
The Kasami sequence set consists of linear combinations of the two 


sequences which, after converting to +1 sequences, are 


, bl+b*. 

a 1 = (-1) t t+ \ 


V t ,i = l,....,2 n/2 -l 


C A6-7) 


and 


?n/2 


= ' H) 
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This yields a set of 



(A6-8) 


sequences, all with period 

L = 2 n -l (A6-9) 

and 

max(C p C 2 ) = 2 n/2 + 1 . (A6-10) 

B. Bent sequences [5,6] posses an underlying arithmetic structure 
in GF(2 n ), n divisible by 1, which is linked to the space V n of binary 
n-tuples over GF(2) by a basis B-j »02»* •• »& n for GF(2 n ) which has the 
property that 

(1 » i = j 

,tr(0.e t ) = (A6-11) 

J 0 s otherwise 


Here tr(-) represents the trace function mapping GF(2 n ) onto GF(2) [7]. 

This generates a correspondence between discrete Fourier transforms 

of functions defined on V and trace transforms of the same functions 

n 

defined on GF(2 n ). This property is exploited along with the fact that 
bent functions [8] on V n have a flat Fourier transform to eventually 
give the following set of sequences: 




(-D 


g(y,)+yTy,+c t x+i t y 
H -1-2 


(A6 -1 2) 


where is the (vector) contents at time t of a Galois-configured 
linear-feedback shift-register with, a primitive characteristic 
polynomial of degree n. 
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where L is a specially designed n/2 x n matrix, the dimension of 
and Y_ 2 ^ ein 9 n /4, I is the representation of i as a binary n/2-tuple 
C^is a fixed non-zero constant and G(_Y^) is a fixed arbitrary Boolean 


Function of . This design results in a set of sequences with the 
same M,L, and max(C-j,C 2 ) parameters as the Kasami sequences (see (8) - 


OO)). 

C. Group character sequences [9,11] are based on properties of the 
group M(L) of integers relatively prime to L under multiplication 
modulo L. In the special case when L is a prime, then 


a 


i 

t 


0 , t = 0 

p U ^, 0 < t < 


where p is a primitive L-l_st root of unity and £(t) 
logarithm of t in the sense that 


( A6-14) 

L, 

is the modulo L 


= t modulo L , (A6-15) 

g being a primitive element of M(L). The index i is restricted to 
1 < i < L-l . Of course for large values of L the use of aj = 1 will 
make little difference in the final results. Group character sequences 
have the following properties: 


H = L-2 

L = prime number 
c i ■ 1 . 

c 2 = 1 /r 

C 11 (0) » L-l 


'(A6-16) 


Generally, the ith sequence is composed of {!_ -1 )/gcd(i ,L-1 ) order roots 
of unity, e.g., i = (L-l )/2 is a sequence of +1 *s and is usually called 
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a quadratic residue sequence. 

D. Welch and Alltop separately have proposed (unpublished to my 
knowledge) signal designs which incorporate a cyclic difference set 
structure [12] to determine the locations of the non-zero elements 
of a sequence. A (v,k,A) cyclic difference set is a collection (t - > 

J 

of integers in the range 0 _< d. < v with the property that the equation 

0 


t. - t. - a mod v 


(A6-17) 


has x solutions for ifO and k solutions when £=0. We view the elements 
of the difference set as the times t at which the sequence elements 
are non-zero. The values of the non-zero elements of each sequence are 
chosen so that any pair of sequences is orthogonal or nearly orthogonal 
when the shift parameter t is zero. 

For example let {x^} be an m-sequence over GF(q), i.e., it satisfies 
an nth order linear recursion over GF(q ) and has period q n -l . Then a 
cyclic difference set with parameters 

v , ail k . aiil 

V q-1 9 K q-1 9 q-1 

is given by 


(A6 -18) 


<& = {t: x t =0> 


(A6-19) 


If we consider a set of sequences based on the above difference set 
with n = 3, then the resulting design parameters are 


-c n (o) = 


L = v = q +q+l 

M = k = q+1 (A6-20) 

max(C-j jCg) = X = 1 . 

Orthogonality of the sequences at t= 0 imposes the result M=k. When q 
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is one less than a multiple of 4, then the rows of a Hadamard 
matrix [13] can supply the appropriate modulation. In this case if 
the L-th entry in the (qfl) x (q+1 ) Hadamard matrix is denoted by b... 

J ‘ 0 

and the elements of the difference set are t-j , . . . jt^, then 

j ( b ij 1 — ’ 1 1 _i J " 1 M 

( 0 otherwise 


(A6-21) 


Other similar designs are possible. 

E. Discrete Linear FM sequences [14,15] of various types have been 
studied. For example. 


J ^ 


- p 


u 2 


0 < t < L, 0 < 1 < p(L) (A6 -22) 


where L is an odd number, p is a primitive L-th_ root of unity, and p(L) 
is the smallest prime divisor of L. When L is prime, 

L = p(L) (A6_23) 


and, using Gaussian sums, it can be shown that 

C = L - ,/2 
max 


(A6-24) 


In addition to this design Alltop [15] suggests optimal designs based 
on phase functions whicn are proportional to t and based on difference 
sets (although sequence elements all have unit magnitude). 

3. Comparisons 

It appears that there are many designs which asymptotically 
achieve the Welch bound on correlation as the sequence period L increases. 
How does one choose a CDMA signal set design from among the class of 
''optimal" designs? 

If you are restricted to binary (+1 ) modulation, the obvious 
candidates are the Kasami sequences and the bent function sequences. 
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The choice may be dictated by L which is 4-1 for Kasami sequences 
k 

and 16 -1 for bent sequences, k being an arbitrary integer in each 
case. The sequences are comparable in terms of implementation com- 
plexity for the same L but the bent sequence set has -two distinct 
disadvantages : 

(1) The Kasami sequences have a linear span on the same order 
of 3n when GF(2 n ) is the basic field, while bent sequences 
apparently have linear spans which can nearly achieve Key's 
upper bound [16], 

Y (?) ( A6 ~ 25 ) 

i=l 


where d is the degree of the bent function (d <_ n/4 
depending on the function chosen). 

(2) All of the Kasami sequences are generated by the same 
hardware, with choice of sequence made by initializing 
register contents. Thus it is difficult to initialize 
a generator to begin producing a copy of a particular Kasami 
sequence at some arbitrary point within the sequence. On 
the other hand bent sequence generators have "time" 
controlled by a shift register and sequence selection 
performed by an independent setting. Hence bent sequence 
generators are easily set to produce a given sequence. 

It is worth noting that Gold codes [17] which are now in use in several 
sys.tems, e.g., [18], have the same drawbacks as Kasami sequences. In 
addition, while Gold codes are a larger -coll ecti on (M=2 n +1) of binary 
sequences for the same period (L=2 n -1), they do not come close to 
achieving the Welch bound. 
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The ability to generate and correlate multiphase sequences consider- 
ably enlarges the variety of periods L for which optimal designs are 
known. The number of distinct phases which must be handled is a function 
of the number of sequences actually required as well as the peri'od 
length. For example the group character sequences of length 1=251 are 
composed in general of 256tji roots of unity, but the ith sequence in the 
set is made up of 256/gcd(256,i )th roots of unity. Hence in this case 
the 128th sequence is the binary quadratic residue sequence, the 64th 

and 192th sequences are composed of 4th^ roots of unity, and in general 
k-1 k 

there are 2 sequences using 2 th_ roots of unity k=l,...,7, with the 
remainder using some primitive 256th roots of unity. 

In comparing the group character sequences with the discrete linear 
FM sequences, one must consider the problem of mechanization for large 
L. The FM sequences have a relatively simple algorithm (A6-22) for 
determining the phase of each bit. On the other hand group character 
sequence generation is based on computing the logarithm of t modulo L-l . 

In most cases this is a difficult computation [19]. 

The main detractions of the difference set design are: (1) irregular- 

ity of the transmitter power, and (2) the requirement of phase coherence, 
despite the on-off nature of the signal. The orthogonal +1 modulation 
of the non-zero pulses is easily achieved, especially when the number k 
of pulses per period is a power of 2. 

Ultimately, with the variety of designs available, it appears that 
convenience of mechanization may very well be the deciding factor in a 
design. 
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APPENDIX 7 

SOLARSIM SUBROUTINE SIDR 

The SOALRSIM subroutine SIDR is developed to compute the Costas 
loop phase jitter as a function of the pilot signal and RF front end 
parameters. Figure A7-1 is a sample run for SIDR. The inputs to 
the program are described as follows: 

F3 - 3 dB bandwidth of the RF filter normalized to the 

chip rate R 
r c 

DELTA - 3 dB bandwidth of the notch filter normalized to R 


GAMMA 
M - 
K1 - 
K2 - 


2 . 


7 is the notch filter attenuation at band center 
PN code length 

Coupling coefficient for noncoherent phase noise 
Coupling coefficient for coherent phase noise 
The other inputs are self explanatory from Fig. A7-1. There are some 
intermediate outputs besides the resultant Costas loop phase jitter 
SIGMA. They are: 

POWER - Received pilot signal power in W 

ALPHA - Time constant of the RF filter normalized to R 

c 

TC - Chip time 

NEQ - Equivalent thermal noise level after despread 

Lp - Fractional loss in pilot signal power 

A software package is also developed to plot the computer output 
(as shown in Figs. 11.1-11.5) on an HP-plotter. 
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ENTER NUMBER OF ELEMENTS OF F3 t BELT A* GAMMA c My RC i 
(ONE FLOATING POINT ENTRY PER LINE) 

1 ♦ 

1 4 

1. 

1 , 

1 * 

ENTER VALUES 

2« 

ENTER VALUES 
♦ I 

ENTER VALUES 

,001 

ENTER VALUES 

10000 , 

ENTER VALUES 

10 , 

ENTER Kt J 

,01 

ENTER l\I? t 

,01 

(FOR THE FOLLOW! NO s’ HIT RETURN FOR NOMINAL VALUES r OTHERWISE ENTER DESIRED VALUE) 
ENTER CHIP DIAS E IN % (NOMINAL 1.0) 

ENTER LOOP BANDWIDTH LN HZ (NOMINAL 10) l 

ENTER TRANSMIT POWER IN KW (NOMINAL 65) i 

ENTER PER CENT OF TRANSMIT EFFICIENCY (NOMINAL 50) i 

ENTER VALUE FOR DIAMETER OF TRANSMIT ANTENNA AREA IN METER (NOMINAL 10) i 

E- 0,1000 BL“ 10,0000 
POWER- 0 , 354033E-06 
kl~ 0,0 1.0000 l\2- 0,010000 

GAMMA = 0.001000 


FOR RF FILTER 3DB BANDWIDTH IN FRACTION OF RC (F3) * 
FOR DELTA l 
FOR GAMMA i 
FUR M i 

FOR CHIP RATE IN MHZ (RC) J 


ALPHAS 0,079577 
DEL T A-' 0 ♦ 9V99996F-01 
M» 10000,0 

RC“ 10. 0000 TC™ 0 ♦ 100000E-06 

NEN« 0 ,441 106E-13 

LF- 0.643650 SIGMA- 0.104654 


Figure A7-1. A Sample Run of the Program SIDR in SQLARSIM. 



